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Interval  analysis  is  applied  to  the  fixed-point  problem  x  *  $(x)  for 
continuous  $ :  s  -*•  S,  where  the  space  S  is  constructed  from  Cartesian  products  of 
the  set  R  of  real  numbers,  with  componentwise  definitions  of  arithmetic 
operations,  ordering,  and  the  product  topology.  Spaces  of  this  kind  include 
many  of  the  usual  Hilbert  and  Banach  spaces  important  in  analysis.  With  the  aid 
of  an  interval  inclusion  *:  IS  ♦  IS  in  the  interval  space  IS  corresponding  to  S, 
the  interval  iteration  process  =  XN  n  ♦  (X^)  is  shown  to  converge  if  the 

initial  interval  XQ  contains  a  fixed  point  x*  of  on  the  other  hand, 

divergence  of  the  iteration  (X^+^  *  0  for  some  N)  proves  that  XQ  contains  no 

fixed  points  of  $,  while  ♦  (Xjj)  c  ^  for  some  N  establishes  the  existence  of  a 

fixed  point  x*  €  and  guarantees  the  convergence  of  the  interval  iteration. 

Each  step  of  interval  iteration  provides  lower  and  upper  bounds  for  fixed  points 
of  $  in  the  initial  interval,  from  which  approximate  values  and  guaranteed  error 
bounds  for  them  can  be  obtained  directly.  In  addition  to  interval  iteration, 
operator  equation  and  dissection  methods  are  also  considered  briefly. 


Since  the  theory  of  interval  iteration  applies  directly  when  only  finite 
subsets  of  S,  IS  are  used,  it  is  adaptable  immediately  to  actual  computation.  A 
numerical  example  is  given  of  the  use  of  interval  iteration  for  the  approximate 
solution  of  a  nonlinear  integral  equation  of  radiative  transfer.  It  is  shown 
that  numerical  results  with  acceptable,  guaranteed  accuracy  can  be  obtained  with 
a  modest  amount  of  computation  for  an  extended  range  of  the  parameter  involved. 


AMS  (MOS)  Subject  Classifications:  47H10,  58C30,  65G10,  65J10,  65J20,  65R20 

Key  words:  Fixed-point  problems,  Interval  iteration,  Error  analysis.  Nonlinear 
integral  equations 


Work  Unit  Number  3  -  Numerical  Analysis  and  Scientific  Computing 


SIGNIFICANCE  AND  EXPLANATION 


Scientific  computation  is  devoted  largely  to  processes  for  numerical 
solution  of  ordinary  and  partial  differential  equations,  integral  equations,  and 
finite  systems  of  linear  and  nonlinear  equations,  all  problems  which  can  be 
formulated  abstractly  as  a  fixed-point  problem  in  some  suitable  space*  Since 
the  computation  is  carried  out  on  a  computer,  and  not  in  the  space  in  which  the 
problem  is  posed,  methods  of  discrete  approximation  are  used,  and  ordinarily 
detailed  analysis  is  required  to  estimate  the  reliability  and  accuracy  of  the 

flood  of  floating-point  numbers  produced,  as  ns  their  _ 

terms  of  the  original  problem.  In  this  report,  which  is  the  text  of  an  invited 
paper  for  the  Special  Session  on  Fixed  Points  and  Operator  Equations  at  the  1983. 
Summer  Meeting  of  the  American  Mathematical  Society  in  Albany,  New  York, 4 it  is 
shown  that  interval  methods  can  be  used  to  overcome  many  of  the  difficulties  of 
this  kind  which  arise  in  scientific  computation.  The  goal  is  to  have  the 
computation  itself  provide  the  required  information  about  the  reliability  and 
accuracy  of  the  computed  results. 

It  is  shown  that  interval  iteration  can  be  used  to  prove  existence  or 
nonexistence  of  solutions  in  given  regions,  and  to  obtain  lower  and  upper  bounds 
for  solutions  in  those  regions.  The  transition  tc  actual  computation  consists 
8 imply  of  identification  of  the  floating-point  numbers  available  on  a  computer 
with  a  finite  set  of  elements  of  the  space  in  question,  for  example,  as  step 
functions  with  floating-point  values,  and  directed  rounding  from  the  space  to 
this  finite  subspace Naturally,  the  computer  used  has  to  support  interval 
arithmetic  by  having  noth  upward  and  downward  rounding  of  arithmetic  operations 
available,  as  required  >y  the  IEEE  Standard  for  floating-point  arithmetic.  This 
capability  is  essential  for  accurate,  reliable  computation. 

A  numerical  example  is  given  in  which  a  nonlinear  integral  equation  of 
radiative  transfer  is  solved  for  an  extended  range  of  parameters  on  a  personal- 
sire  microcomputer  in  only  a  few  seconds.  For  suitable  initial  intervals,  this 
modest  computation  gives  a  guarantee  of  the  existence  of  the  solution  of  the 
integral  equation,  and  proves  that  the  numerical  results  obtained  are  accurate 
to  better  than  10%.  In  other  cases,  the  computation  shows  that  the  initial 
interval  is  not  large  enough  to  contain  a  solution  of  the  integral  equation, 
which  is  also  valuable  information  when  exploring  for  regions  in  which  solutions 
will  lie. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC ,  and  not  with  the  author  of  this  report. 


INTERVAL  METHODS  FOR  FIXED-POINT  PROBLEMS 

L.  B.  Rail 

1 .  Real  spaces .  Here,  the  fixed-point  problem 

(1.1)  x  ■  $(x) 

will  be  considered  for  continuous  operators  S  ♦  S  acting  In  spaces  S  which 
can  be  constructed  In  a  simple  way  from  the  set  R  of  real  numbers.  This 
construction  is  based  on  the  Cartesian  product 

(1.2)  P  -  J[  R  , 

of  the  real  numbers  over  an  Index  set  A.  Elements  f  of  P  are  vectors  with  real 
ccMponents  denoted  by  fa  or  f(a),  a  e  A. 

For  example,  if  A  ■  {1,2,...,n}  is  the  set  of  the  first  n  positive 
integers,  then  P  is  the  set  Rn  of  n-dimensional  real  vectors  f  -  (f  ,,,f  2, . . . fR) » 
if  A  -  fa,b),  an  interval,  then  P  is  the  set  R[a,b]  of  all  real-valued  functions 
f  on  [a,b],  with  fa  ■  f(a). 

In  general,  the  product  set  P  can  thus  be  considered  to  be  the  set  of  all 

real-valued  functions  on  A,  P  »  (f  |  fs  A  ♦  R)  [1]. 

\ 

The  set  P  will  now  be  equipped  with  arithmetic  operations,  a  partial 
ordering,  and  a  topology,  all  taken  from  the  underlying  set  R  of  real  numbers. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041. 


Definition  1.1.  The  product  set  P  given  by  (1.1)  is  called  a  basic  real 
space  over  A  if: 

(i)  The  arithmetic  operations  +  ,  -  ,  •  ,  /  are  defined  in  P  by 

M.3)  (f  *  g)a  =  fQ  *  go,  a  e  A, 

for  *  e  {+,  -,  *,  /} ,  with  divisors  restricted  to  the  set 

H.4)  D  =  {g  |  gfl  ¥  0,  a  e  A>: 

(ii)  P  is  partially  ordered  by  _<_,  where 

(1.5)  f  <_  g  iff  fa  _<_  ga,  a  e  A; 

(iii)  The  topology  of  P  is  the  product  topology  of  the  ordinary 
(metric)  topology  of  the  reals  [12]. 

Properties  of  arithmetic  operations  in  a  basic  real  space  P,  such  as 
commutivity  of  addition  and  multiplication,  thus  follow  from  the  corresponding 
properties  of  real  arithmetic,  and  the  partial  ordering  induced  in  P  by  the 
order  relationship  for  real  numbers  results  in  P  being  a  complete  lattice  [5] . 

A  basic  real  space  will  be  a  linear  space  if  multiplication  of  a  vector 
f  «  P  by  a  scalar  a  €  R  is  defined  by  multiplication  of  f  by  the  vector  in  P 
with  all  components  equal  to  a.  Vectors  which  have  all  components  equal  (that 
is,  which  are  constant  on  A),  will  be  denoted  on  occasion  by  this  common  value, 
in  particular,  the  notation  0  and  1  will  be  used  for  the  vectors  such  that  * 

0  and  1ffl  •  1  for  all  a  £  A,  respectively.  Thus,  for  example,  for  A  =  [0,1],  the 
basic  real  space  P  3  R[0,1]  includes  the  elements  of  the  Banach  space  C[0,1]  of 
continuous  real  functions  on  [0,1],  and  the  Hilbert  space  consisting  of 
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all  real  functions  which  are  square- integr able  in  the  sense  of  Lebesgue  on 

[0,1]. 


Instead  of  the  entire  space  P,  attention  can  be  restricted  to  the  subspace 
of  P  consisting  of  bounded  elements  f,  for  which 

(1.6)  Ifl  -  aup  |f  |  <  +  •  . 

a  £  A  a 

On  this  subspace,  divisors  have  to  be  restricted  to  the  set  of  elements  g  which 
have  bounded  reciprocals  1/g. 

The  absolute  value  |  |  and  other  functions  of  elements  of  basic  real  spaces 
are  again  defined  componentwise,  provided  this  makes  sense  for  the  given 
function. 

In  order  to  consider  higher-dimensional  problems,  a  more  general  class  of 
spaces  is  introduced  by  the  following  recursive  definition. 

Definition  1.2.  A  zeal  space  S  is  either  a  basic  real  space  or  the 
Cartesian  product  of  real  spaces.  In  the  latter  case,  arithmetic  operations  and 
partial  ordering  of  the  product  space  are  defined  componentwise  on  the  basis  of 
the  corresponding  operations  and  relations  in  the  factor  spaces,  and  the 
topology  for  the  product  space  is  the  product  of  the  topologies  of  its  factors. 

For  example ,  the  real  space  S  ■  Rn  x  Rm  consists  of  all  functions  which  map 
the  n-dimensional  space  Rn  into  the  m-dimensional  space  Rra. 

Real  spaces  will  sometimes  be  called  R-spacee  for  brevity. 

As  is  generally  known,  an  overwhelming  number  of  theoretical  and  applied 
mathematical  problems  can  be  expressed  in  the  deceptively  simple  form  (1.1).  In 
this  paper,  this  general  problem  will  be  discussed  from  the  standpoint  of 
interval  analysis  111},  [18],  [3],  [4].  The  next  sections  will  introduce  the 
necessary  ideas  of  interval  spaces  and  operators. 
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2.  interval  spaces.  Given  a  real  space  S,  the  corresponding  interval 


space  IS  is  constructed  in  the  same  way  as  the  set  IR  of  nonenpty,  closed, 
bounded  intervals 

(2.1)  [a,b]  *{x  |  -•<a<xi<b<  +•} 

is  obtained  from  the  set  of  real  numbers  R.  In  IR,  the  operations  of  interval 
aritbaetic  [17],  [18]  are  defined  by 

(2.2)  [a,b]  *  [c,d]  -{z|z»x*y,  xe  [a,b] ,  y  c  [c,d]}, 

*  e  {+,  -,  *,  /},  with  divisors  restricted  to  the  set 

(2.3)  D  -  { [c,d]  |  cd  >  0), 

that  is,  to  the  set  of  intervals  which  do  not  contain  zero.  Intervals  [c,d] 
belonging  to  D  will  have  bounded  reciprocals 

(2.4)  [c,d]-1  -  [1/d  ,  1/c]. 


It  is  important  to  note  that  intervals  do  not  form  a  group  with  respect  to 
addition,  since  additive  inverses  do  not  exist  in  general;  for  example,  one  has 
[0,1]  -  [0,1]  -  [-1,1]  instead  of  [0,0]. 

A  partial  ordering  can  be  defined  in  IR  by,  for  example 

(2.5)  [a,b]  ±  [c,d]  iff  b  <_  c, 
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([17],  p.  7).  There  is  also  a  metric  topology  for  intervals  [17].  Use  of  these 
ideas  leads  to  the  following  definition. 


Definition  2.1.  The  interval  apace  IS  corresponding  to  a  real  space  S  is 
the  set  of  intervals  [a,b]  analogous  to  (2.1),  where  a,b  are  bounded  elements  of 
S  and  a  £  b  In  the  partial  ordering  on  S,  with  componentwise  definition  of 
interval  arithmetic  and  partial  ordering  based  on  (2.2)  and  (2.5)  respectively, 
and  topologized  by  the  product  of  the  metric  topology  for  intervals  ([17],  Chap. 
4). 

In  IS,  division  is  restricted  to  intervals  with  bounded  reciprocals* 
Elements  x  of  S  can  be  identified  with  the  degenerate  intervals  [x,x]  e  IS 
having  equal  lower  and  upper  endpoints.  The  notation 

(2.6)  x  *  [x,x] 

will  be  used  throughout  to  indicate  this  identification. 

The  interval  space  IS  corresponding  to  a  real  space  S  can  also  be 
constructed  from  the  intervals  IR  of  real  numbers  R  according  to  the  recipe  for 
S,  keeping  only  intervals  with  bounded  endpoints.  Although  the  real  space  S  is 
a  linear  space,  the  corresponding  interval  space  IS  will  not  be  a  linear  space, 
because  interval  addition  does  not  have  the  group  property,  as  noted  above. 

Since  intervals  are  sets,  the  inclusion  relation  c  is  defined  in  IR,  and 
provides  another  partial  ordering  of  this  set  [5].  Furthermore,  the 
intersection  [a,b]  n  [c,d]  of  two  intervals  in  IR  is  either  an  interval  or  the 
empty  set  0.  These  ideas  extend  componentwise  to  interval  spaces. 

Definition  2.2.  If  F,G  €  IS,  where  IS  is  an  interval  space,  then  F  c  g 

means 

(2.7)  F  c  G  ,  a  e  A. 

a  a 
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Definition  2.3.  If  F,G  f  IS,  an  interval  space,  and  FQ  n  Gffl  *  ft  for  some 

a  e  A,  then  the  intersection  of  F  and  G  is  said  to  be  empty,  and  one  writes 

(2.8)  F  n  G  -  0. 

Thus,  the  intersection  of  elements  of  an  interval  space  will  either  be  an 

element  of  that  space  or  empty.  On  occasion,  elements  of  an  interval  space  IS 

will  be  referred  to  simply  as  intervals. 

Since  the  intersection  of  a  nested  sequence  of  closed,  nonempty  subsets  of 
R  is  nonempty  [33],  it  follows  that  for  intervals  X^  e  ir  such  that 

(2.9)  x,  a  X2  ^  X3  =  ...  a  Xn  a  ..., 

one  has 

Oft 

(2.10)  n  X  f  0, 

„  n 
n*»0 

since  elements  of  IR  are  closed  and  nonempty  by  definition.  It  follows  from 
Definition  2.4  by  construction  and  the  Cantor  theorem  [33]  that  the  same  result 
holds  in  interval  spaces  IS,  that  is,  if  (2.9)  holds  for  X^  £  IS,  i*  1,2,3,..., 
then  (2.10)  is  true  [28]. 

3.  Interval  inclusions  of  real  operators.  Suppose  that  S,T  are  real 
spaces,  and  IS,  IT  are  the  corresponding  interval  spaces.  An  operator  f:  S  ♦  T 
will  be  called  a  real  operator,  and  an  operator  F:  IS  ♦  IT  will  be  called  an 
interval  operator.  Although  real  operators  can  be  classified  as  linear  or 
nonlinear,  this  cannot  be  done  for  interval  operators,  since  the  interval  spaces 
IS,  IT  are  not  linear  spaces,  and  the  definition  of  linearity  for  an  operator 
depends  on  a  linear  substructure  for  the  spaces  involved. 
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Definition  3.1.  The  interval  operator  F:  IS  ♦  IT  is  said  to  be  ( inclusion) 


monotone  if  for  X,Z  €  is, 

(3.1)  X  c  z  implies  F(X)  c  F(z). 

Definition  3.2.  The  interval  operator  F:  IS  *  IT  Lb  said  to  oontain  the 
real  operator  f:  S  ♦  T  if 

(3.2)  f(X)  =  {y  |  y  =  f(x),  x  e  X}  c  f(X),  X  £  IS. 

This  is  symbolized  by  f  £  F,  and  F  is  said  to  be  an  interval  inclusion  of 
f.  If  F  is  monotone,  then  it  is  called  a  monotone  inclusion  of  f. 

Definition  3.3.  The  real  operator  f:  S  +  T  is  said  to  be  a  restriction  of 
the  interval  operator  F:  IS  ♦  IT  if 

(3.3)  f(x)  =  F(x)  =  F([x,x]),  x  e  X, 

in  the  sense  of  the  identification  (2.6)  of  elements  of  S  with  degenerate 
intervals  in  IS.  If  (3.3)  holds  for  a  (monotone)  inclusion  F  of  f,  then  F  is 
called  a  (monotone)  interval  extension  of  f  [17]. 

Interval  extensions  of  rational  functions  can  be  constructed  by  the  use  of 
interval  arithmetic  [17],  [18],  13],  [4].  The  class  of  operators  for  which 
interval  extensions  are  available  can  be  broadened  by  the  use  of  interval 
extensions  of  elementary  functions,  and  so  on.  Integral  operators  can  be 
extended  by  the  use  of  interval  integration  [9] ,  [26] .  In  actual  numerical 
computation,  it  is  usually  impossible  to  realize  the  restriction  property  (3.3); 
however,  interval  arithmetic  with  directed  rounding  can  be  used  to  construct 
interval  inclusions  of  real  operators  [15],  [17],  [18].  The  results  given  below 


require  only  the  availability  of  interval  inclusions  F  of  real  operators  f,  and 
do  not  depend  on  the  restriction  property. 

4.  Interval  iteration.  Suppose  that  S  *  S  is  a  continuous  operator  in 
a  real  space  S,  and  4:  IS  ♦  IS  is  an  interval  inclusion  of  4. 

Definition  4.1.  The  sequence  {Xn)  of  elements  of  IS  defined  by 

(4.1)  Xn+1  =  Xn  °  *(Xn)'  "=0,1,2,..., 

is  said  to  be  generated  by  interval  iteration  with  initial  interval  XQ  and 
iteration  operator  4. 

Definition  4.2.  If  a  positive  integer  N  exists  such  that 

(4.2)  n  *0^)  =  0, 

then  the  interval  iteration  (4.1)  is  said  to  diverge;  otherwise,  it  converges  to 
the  limit 

09 

(4.3)  X*  =  n  X  €  IP. 

n=0  ** 

The  existence  of  X*  in  the  case  of  convergence  follows  from  the  theorem  of 
Cantor  [33]  cited  in  §2,  since  the  iteration  sequence  (xn)  is  a  nested  sequence 
(2.9)  of  nonempty  intervals  in  S.  If  the  iteration  diverges,  then  the  iteration 
sequence  will  consist  only  of  the  N  intervals  Xq,  X^,  . ..,  X^.  The  following 
theorems  relate  the  convergence  or  divergence  of  the  interval  iteration  (4.1)  to 
the  existence  or  nonexistence  of  a  fixed  point  of  4  I"  the  initial  interval  Xq 
[28]. 
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Theorem  4.1.  If  the  initial  interval  Xq  contains  a  fixed  point  of  $,  then 
the  interval  iteration  (4.1)  will  converge  to  a  limit  X*  which  contains  all 
fixed  points  x*  of  $  which  belong  Xq. 

The  proof  of  this  theorem  [28]  follows  from  the  simple  observation  that  if 
x*  £  Xn  is  a  fixed  point  of  $,  then  x*  £  *(Xn),  so  **  e  Vi'  and  thus ,  by 
mathematical  induction,  the  elements  of  the  iteration  sequence  {Xn }  will  be 
nonempty  and  contain  all  fixed  points  of  $  which  belong  to  XQ  . 

Theorem  4.2.  If  the  interval  iteration  (4.1)  diverges,  then  *:he  initial 
interval  XQ  does  not  contain  a  fixed  point  x*  of  $. 

This  is  the  contrapositive  of  Theorem  4.1. 

Divergence  of  interval  iteration  will  be  observed  in  a  finite  number  of 
steps  if  it  occurs.  On  the  other  hand,  interval  iteration  will  converge  if  XQ 
is  large  enough  to  include  at  least  one  fixed  point  of  $.  Interval  and  ordinary 
iteration  will  be  con^jared  in  more  detail  in  the  next  section.  The  next 
theorems  give  sufficient  conditions  for  the  convergence  of  interval  iteration. 

Theorem  4.3.  If 

(4.4)  *(XN)  c  Xjj 

for  some  positive  integer  N,  then  a  fixed  point  x*  e  X^  c  xQ  of  $  exists,  and 
the  interval  iteration  (5.1)  converges  to  the  limit  X*  given  by  (4.3),  and 
x*  e  X*. 

The  existence  of  x* e  Xj^  follows  from  (4.4)  by  the  Schauder  fixed  point 
theorem  [32],  since  $  is  assumed  to  be  continuous,  and,  as  a  subset  of  S,  X^  is 
closed,  convex  [29] ,  and  compact.  The  convergence  of  the  the  interval  iteration 
to  the  limit  X*  with  x*  <r  X*  then  follows  from  Theorem  4.2. 


Theorem  4.4.  If 


(4.5)  Xj,  c  *0^) 

for  some  positive  integer  N,  then  the  interval  iteration  (4.1)  converges  to  the 
limit  X*  = 

It  follows  from  (4.1)  and  (4.5)  that  X^^  =  X^,  and  thus  Xn  »  X^  for  n  >_  N, 
which  gives  X*  =  by  (4.3). 

Theorems  4.3  and  4.4  show  that  it  is  sufficient  for  the  convergence  of 
interval  iteration  for  X^  and  4  (X^ )  to  be  comparable  for  the  inclusion 
relation  <=  for  some  positive  integer  N. 

Note  that  convergence  of  the  interval  iteration  under  condition  (4.5)  does 
not  imply  the  existence  of  a  fixed  point  x*  e  Xj^  of  4.  It  is  easy  to  construct 
examples  in  which  (4.5)  holds,  but  the  initial  interval  XQ  does  not  contain  a 
fixed  point  of  4*  In  this  case,  it  seems  more  appropriate  to  say  that  the 
interval  iteration  has  "stalled"  instead  of  "converged". 

An  example  of  the  application  of  interval  iteration  to  a  nonlinear  integral 
equation  will  be  given  below. 

5.  Comparison  with  ordinary  iteration.  A  customary  approach  to  the  fixed 
point  pro'*1 -am  (1.1)  is  by  ordinary  iteration,  that  is,  generation  of  the 
iteration  sequence  {x.j}  by 

(5.1)  xn+1  =  <Mxn>»  n  =  0,1,2,..., 

starting  from  the  initial  point  Xq. 

It  is  well-known  that  if  4  is  continuous  and  the  iteration  sequence 
converges  to  x*,  then  x*  is  a  fixed  point  of  4>  conversely,  if  4  does  not  have  a 
fixed  point,  then  the  iteration  sequence  is  divergent  in  the  sense  that  it 
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cannot  converge  to  an  element  of  the  space.  Comparison  of  these  results  with 
those  of  §4  for  interval  iteration  indicc.ts  that  the  conditions  relating 
existence  or  nonexistence  of  fixed  points  to  convergence  or  divergence  of  the 
corresponding  iteration  sequences  are  converse  in  a  certain  sense  [28] .  For 
example,  existence  of  a  fixed  point  x*  of  ^  in  the  initial  interval  Xg  is  a 
sufficient  condition  for  the  convergence  of  interval  iteration  (Theorem  4.1), 
while  this  is  a  necessary  condition  for  the  convergence  of  ordinary  iteration. 
Similarly,  divergence  of  ordinary  iteration  is  a  necessary  condition  for  the 
nonexistence  of  a  fixed  point,  while  divergence  of  interval  iteration  is  a 
sufficient  condition  for  the  nonexistence  of  x*  in  the  initial  interval  XQ . 

Thus,  interval  and  ordinary  iteration  complement  each  other  in  a  certain 
sense,  and  one  or  the  other  can  be  used  as  appropriate  to  a  given  fixed  point 
problem.  One  advantage  of  interval  iteration  for  numerical  purposes  is  that  it 
provides  lower  and  upper  bounds  for  fixed  points  at  each  step.  One  can  write 

(5.1)  Xn  =  [an,  bn],  n  =  0,  1,  2,  ..., 
where  aR,  bn  €  S.  For  each  fixed  point  x*  e  Xg  of  $,  one  has 

(5.2)  ao  .£  ai  <.  a„  <.  •  •  •  <.**.£*••  1.  bn  —  *  *  ’  —  b1  —  b0 , 

and  in  the  limit,  for  X*  =*  [a*,  b*] , 

(5.3)  a*  <  x*  <.  b*. 

Thus,  if  one  knows  initially  that  Xg  contains  at  least  one  fixed  point  x*  of  $, 
then  interval  iteration  provides  a  way  to  obtain  improved  approximations  to  x*, 
with  guaranteed  componentwise  error  bounds.  For 


(5.4) 


where  the  norm  is  defined  by  (1.6). 

Error  bounds  for  ordinary  iteration  (5.1),  on  the  other  hand,  are  not 
usually  obtained  from  the  iteration  itself,  but  require  the  construction  of  some 
majorant  function  E(n)  such  that 


(5.7)  lx*  -  x  I  <  E(n),  n  =  0,1,2,...  . 

n  ~ 

Another  distinction  between  ordinary  and  interval  iteration  arises  in 

actual  computation,  due  to  the  fact  that  the  sequence  {x^}  cannot  be  obtained 

exactly,  but  has  to  be  approximated  by  a  sequence  {z^} .  For  ordinary  iteration, 

this  means  that  lx*  -  z  I  has  to  be  estimated  as  a  separate  issue,  for  example, 

n 

by  use  of  (5.7)  and  a  bound  for  the  approximation  error  lx  -  z  I.  Thus,  the 

n  n 

additional  analysis  which  has  to  be  brought  to  bear  on  the  error  estimation 
problem  can  be  extensive.  In  the  case  of  interval  iteration,  however,  the  lower 
and  upper  bounds  («n)»  {bn}  are  simply  selected  from  a  finite  subset  of  of  S, 
and  the  results  (5.2)-(5.6)  hold  without  modification.  Furthermore,  this  means 
that  a  positive  integer  N  will  exist  for  which  one  observes  either  (4.2) 
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(divergence)  or  (4.5)  (convergence).  Hence,  convergence  or  divergence  of 
Interval  Iteration  Is  observable  in  practice,  at  least  in  principle,  while  the 
approximate  iteration  sequence  (zn)  obtained  by  ordinary  iteration  can  wander 
indefinitely  without  revealing  whether  it  is  convergent  or  not. 

It  is  also  usual  in  application  of  ordinary  iteration  that  the  initial 
point  Xg  has  to  be  chosen  "close  enough"  to  x*  to  obtain  convergence,  and, 
indeed,  some  fixed  points  repel  the  iteration  sequence  {x^}  for  Xg  arbitrarily 
close  to  x*.  As  mentioned  above,  the  initial  interval  Xg  simply  has  to  be  "big 
enough"  to  include  x*  for  interval  iteration  to  converge.  Of  course,  the 
interval  iteration  may  not  converge  to  anything  useful  if  its  limit  X*  is  too 
large. 

The  challenge  in  interval  iteration  is  to  construct  the  interval  inclusion 
•  of  4  in  such  a  way  that  4 (X )  is  not  grossly  larger  than  4(X),  so  that 
approximations  and  error  bounds  which  are  obtained  from  the  computation  itself 
will  be  as  accurate  as  possible.  In  particular,  if  $  is  a  contraction  operator 
in  the  ordinary  sense,  one  would  want  4  to  be  an  Interval  contraction  [6] ,  or  as 
close  to  one  as  can  be  obtained  with  outward  rounding.  Mean  value  and  Taylor 
forms  [7] ,  [29]  have  been  found  to  be  very  useful  for  the  construction  of 
accurate  interval  inclusions  of  sufficiently  differentiable  operators  in  real 
spaces . 

Before  leaving  the  subject  of  interval  iteration,  it  should  be  noted  that  a 
substantial  number  of  the  1200  references  in  a  recent  bibliography  on  interval 
mathematics  [12]  deal  with  some  aspects  of  the  theory  or  examples  of  application 
of  this  technique.  Most  of  these  studies  focus  on  conditions  under  which  the 
inclusion  (4.4)  will  hold,  so  that  the  interval  iteration  will  yield  a  guarantee 
of  existence  of  the  fixed  point,  as  well  as  lower  and  upper  bounds  for  it. 
However,  as  pointed  out  above,  interval  iteration  can  also  be  useful  to  prove 
nonexistence  of  fixed  points,  as  well  as  to  obtain  bounds  for  fixed  points  known 
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to  belong  to  Xq  on  the  basis  of  some  other  argument.  In  addition  to  the 
pioneering  work  of  Moore  [17],  [18]  substantial  contributions  by  Alefeld  [2], 
Krawczyk  [13],  [14],  Nickel  [20],  [21],  and  wisskirchen  [34]  should  be  noted, 
among  others.  The  books  by  Alefeld  and  Herzberger  [3] ,  [4]  also  contain  much 
valuable  information  on  this  topic.  The  interval  bibliography  [11]  also 
contains  references  to  proceedings  of  conferences  on  recent  developments  in  this 
rapidly  growing  field. 

6.  Operator  equation  methods.  Interval  methods  can  also  be  based  on  the 
equivalence  of  the  fixed  point  problem  (1.1)  to  an  operator  equation 

(6.1)  f(x)  -  0, 

where  f:  S  ♦  S  is  such  that  f(x*)  *  0  if  and  only  if  x*  *  4(x*).  Such 
equivalent  problems  can  be  constructed  in  many  ways,  for  example,  simply  by 
taking  f  to  be  defined  fcy 

(6.2)  f(x)  =  x  -  $(x). 

In  many  cases,  an  equivalence  between  (1.1)  and  (6.1)  is  exploited  in  the 
other  direction;  one  attempts  to  solve  (6.1)  for  x  ■  x*  by  constructing  an 
iteration  operator  4  for  which  x*  is  a  fixed  point,  and  then  using  ordinary 
iteration  (5.1).  Typical  examples  of  this  are  furnished  by  the  use  of  Newton's 
methods  and  its  variants  for  the  solution  of  operator  equations  [24] . 

Suppose  that  Fs  IS  ♦  IS  is  an  interval  inclusion  of  the  operator  f:  S  ♦  S 
appearing  in  (6.1).  In  searching  for  an  initial  interval  Xq  which  contains  a 
fixed  point  x*  of  an  equivalent  operator  4/  the  following  exclusion  theorem  is 
often  useful. 
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Theorem  6.1.  If 


(6.3)  0  i  F(X), 

then  the  interval  X  does  not  contain  a  solution  x  ■  x*  of  equation  (6.1). 

Since  F  is  an  inclusion  of  f,  (6.3)  contradicts  the  existence  of  an  element 
x*  6  X  such  that  f(x*)  ■  0. 

7.  Dissection  methods.  The  search  for  fixed  points  of  the  operator 
S  ♦  S  in  intervals  X  e  IS  can  be  carried  out  with  the  aid  of  a  three- valued 
logic  A  »  {TRUE,  7,  FALSE}  and  a  test  T^  such  that 

(a)  T,  (X)  “  TRUE  means  X  contains  a  fixed  point  x*  of  9> 

9 

(7.1)  (b)  T . (X )  ■  7  means  the  test  is  inconclusive. 

9 

(c)  T^(X)  «  FALSE  means  X  does  not  contain  a  fixed  point  x*  of 

As  shown  above  in  §4,  interval  iteration  starting  from  XQ  =  X  provides  an 
example  of  such  a  test,  with  (7.1) (a)  holding  if  (4.4)  is  observed,  (7.1 )(c)  if 

(4.2)  is  observed,  and  no  conclusion  about  existence  or  nonexistence  of  x*  can 
be  drawn  if  the  interval  iteration  simply  converges  without  (4.4)  being  true  at 
some  step,  for  example,  if  (4.5)  occurs  with  proper  inclusion.  In  actual 
computation  on  finite  subsets  Sg  and  ISg  of  S  and  IS,  one  of  these  alternatives 
will  be  observed  in  a  finite  number  of  steps,  so  that  interval  iteration  can  be 
used  as  the  test  T^.  Of  course,  other  criteria  for  existence  or  nonexistence  of 
fixed  points  x*  of  $  in  X  can  be  applied  instead  of  interval  iteration,  or  to 
supplement  it. 

Suppose  now  that  X  is  an  interval  such  that  T^(X)  -  7,  and  it  is  desired  to 
determine  whether  or  not  X  actually  contains  a  fixed  point  x*  of  $.  Then,  X  can 
be  dissected  into  intervals  X^,  X*2*,  x^,  such  that 
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(7.2) 


X 


r(  1  ) 


U  x 


(2) 


t(k>. 


(The  usual  case  k  *  2  is  called  bisection  of  X. )  The  test  T^  can  now  be  applied 
to  the  X*1-*,  i  -  1,2,...,k,  in  turn.  If  no  positive  conclusion  is  reached,  then 
the  subintervals  X^  for  which  T^(x^)  “  FALSE  can  be  discarded,  and 
dissection  applied  again  to  the  remaining  intervals  until  one  is  found  which 
contains  a  fixed  point  of  or  all  are  rejected.  In  the  latter  case,  of 

course,  the  original  interval  X  is  fixed-point  free. 

Once  a  subinterval  has  been  found  which  contains  fixed  points,  then 
interval  iteration  or  some  other  method  can  be  used  with  confidence  to  obtain 
their  values,  at  least  approximately.  Dissection  is  most  appropriate  in  finite- 
dimensional  spaces,  however,  even  in  spaces  of  low  dimension,  considerable 
bookkeeping  can  be  required  to  keep  track  of  the  subintervals  still  under 
consideration,  and  is  thus  a  task  suitable  for  a  computer.  Furthermore,  it  has 
been  shown  by  Hoore  and  Jones  [19]  that  bisection  will  yield  a  definite 
conclusion  (7.1) (a)  or  (7.1)(c)  for  certain  fixed  point  problems  in  a  finite 
number  of  steps.  In  particular,  there  are  classes  of  fixed  points  for  which 
well-known  sufficient  conditions  for  their  existence  hold  in  regions  surrounding 
the  fixed  point,  and  these  regions  contain  intervals  in  which  the  fixed  point 
lies.  In  this  case,  a  positive  conclusion  will  be  reached  after  the  dissection 
has  been  carried  out  to  the  point  that  the  subintervals  under  consideration  are 
sufficiently  small.  An  example  of  this  situation  is  furnished  by  the  famous 
theorem  of  Kantorovich  on  the  convergence  of  Newton's  method  [24],  [23],  [25]. 

8.  A  nonlinear  integral  equation.  The  above  results  will  be  illustrated 
by  application  to  the  numerical  solution  of  the  integral  equation 


(8.1) 


H(m> 


1 _ 

X  ,1  uH(u' ) 
2  j  U  +  M' 


du' 


0  <  p  <  1, 
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which  is  a  fora  of  tha  H-equation  of  radiative  transfer  [10].  Equation  (8.1) 
will  be  considered  to  be  the  fixed  point  problem 


(8.2)  H  -  $ (H ) 

in  the  space  R[0 , 1]  of  real-valued  functions  on  the  interval  [0,1].  For 
numerical  purposes,  the  subset  R^[0,1]  will  be  taken  to  the  set  of  step 
functions  with  constant  floating-point  values  on  the  subintervals 


(8.3) 


-  <i-1)/10  ini.  i/10  -  vL,  i  -  1,2,..., 10. 


This  space  and  the  corresponding  interval  space  IR^ [0,1]  give  a  "low  resolution" 
view  of  R[0, 1 ]  and  IR[0,1],  respectively .  For  confutations 1  purposes,  these 
spaces  of  piecewise  constant  functions  and  intervals  can  be  identified  with  the 
vector  space  R#10  of  vectors  x  ”  (x1 ,x2, . . . ,x10 )  with  floating-point  numbers  as 
components ,  and  the  corresponding  interval  space  XR^10  of  interval  vectors 
X  ■  (X1 ,X2, . . . ,X1Q ),  respectively. 

An  interval  inclusion  9  of  the  operator  $  defined  by  the  right  side  of 
(8.1)  will  be  constructed  with  the  aid  of  interval  integration  [9],  [26],  as  in 
[31].  Identifying  an  element  x  of  R#[0,1]  by  the  vector  x  in  R#10  of  its 
values,  the  integral  transform  of  x  is  given  by 


(8.4) 


I 

0 


1  UxlUll 

y  +  y' 


dy ' 


1° 

S  x,ln{(y  +  y.)/(y  +  y.  )}, 
j»1  3  3  3 


yi_1  <.  y  <.  y^,  i  •  1,2,.. .,10.  in  each  subinterval  , u^3 »  the  function 

(8.4)  is  rounded  downward  or  upward  to  a  floating-point  constant  to  obtain  an 
element  I(x)  of  R^[0,1].  Applied  to  the  lower  and  upper  endpoint  functions  of 
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an  element  X  of  1R|[0,1],  this  technique  yields  an  inclusion  I(X)  e  IR$[0,1]  of 
the  interval  integral  transform  of  an  element  X  £  IR#  [0,1].  The  rest  of  the 
construction  of  4:  IR#[0,1]  ♦  IR#[0,1]  defined  by 

(8*5)  *(X)  ■ - 1 -  ,  x  £  XR„  [0,1], 

1  -  f  I(X)  * 

is  now  done  by  rounded  interval  arithmetic  in  R^10  [4],  [15],  [17],  [18],  with 
the  division  in  (8.5)  carried  out  componentwise. 

The  range  of  values  of  X  of  physical  interest  is  0  <_  X  £  1  [10].  By 
ordinary  methods  of  functional  analysis,  it  is  easy  to  establish  that  equation 
(8.1)  has  a  solution  for  0  X  <_  1/(2  In  2)  <•  0.721...  [24].  The  same  bound  has 
been  obtained  for  interval  iteration  applied  to  the  standard  form  of  (8.1), 
using  a  cruder  form  of  interval  integration  than  given  here  [22] ,  the  basic 
result  being  that  Theorem  4.3  holds  for  the  initial  interval 

(8.6)  X  -  [  1  ,  -  -~n~ -  {1  -  /I  -  2  X  In  2  }  ]. 

Here,  interval  iteration  will  be  applied  to  (8.1)  for  0.74  <.  X  <_  1.00  in 
order  to  extend  the  range  of  X  beyond  that  established  previously.  The 
iteration  is  performed  in  the  "Gauss-Seidel"  mode  recommended  by  Caprani  and 
Madsen  [8]  for  integral  equations,  in  which  as  much  updated  information  is  used 
at  each  step  as  possible,  that  is,  the  _ith  component  $(X)^  of  4(X)  is  computed 
by 

(8.7)  KX^  -  *(*(X)1,...,*(X)i_1,Xi,Xi+1,...,X10), 

i  -  1,2,...,  10.  Wisskirchen  [34]  has  established  the  significant  result  that 
interval  iteration  performed  in  this  way  converges  to  the  same  limit  as  when 
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straightforward  evaluation  of  ♦(X1 ,X2, . • • jX^ )  is  used.  However,  since  (8.7) 
uses  smaller  component  intervals  at  each  substage  of  the  evaluation  of  4(X),  one 
can  expect  more  rapid  convergence.  (This  is  another  distinction  between 
interval  and  ordinary  iteration,  since  a  Gauss-Seidel  version  of  a  convergent 
ordinary  iteration  method  may  diverge  or  converge  to  a  different  limit. ) 

The  results  are  given  in  Table  8.1.  For  each  initial  interval  and  values 
of  X  listed,  the  three  following  columns  report  the  number  N  +  1  of  the 
iteration  at  which  (4.4)  is  observed  (Existence),  successive  intervals  agree  to 
five  decimal  places  (Convergence),  or  (4.2)  happens  (Nonexistence).  Convergence 
in  the  sense  of  Theorem  4.4  will,  of  course,  follow  (4.4)  after  a  certain  number 
of  iterations,  but  can  also  occur  without  a  guarantee  of  existence  of  a  solution 
if  the  inclusion  in  (4.S)  is  proper. 

Examination  of  the  given  results  reveals  that  relatively  few  iterations  are 
required  to  establish  existence  or  nonexistence  of  a  fixed  point  in  the  initial 
interval  when  one  of  these  outcomes  obtains.  The  number  of  iterations  needed 
for  convergence,  on  the  other  hand,  depends  on  the  number  of  digits  to  which  the 
endpoints  of  successive  intervals  are  required  to  agree,  that  is,  on  the 
precision  of  the  floating-point  number  system  used.  Experiments  with  this 
example  indicate  that  this  dependence  is  linear,  with  twice  as  many  iterations 
needed  to  obtain  convergence  to  twelve  significant  decimal  digits  as  were 
observed  for  six.  Since  step  functions  are  relatively  crude  approximations  to 
continuous  functions,  particularly  on  the  coarse  subdivision  of  [0,1]  into  only 
10  subintervals,  it  is  a  waste  of  time  to  insist  that  this  interval  iteration 
converge  to  more  than  a  few  decimal  places. 

As  observe!  in  [22],  [31],  better  approximations  to  the  solution  of  the 
integral  equation  are  obtained  at  the  break  points  ■  (0.1 )i,  i  -  1,2,..., 9, 
than  in  the  corresponding  subintervals,  since  the  solution  evaluated  at  these 
points  has  to  satisfy 
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1 


(8.8) 


x*<»±>  £  *!  n  Xi+1  -  Zit  i  -  1,2 . 9. 

An  approximation  to  x*(y^)  can  be  taken  to  be  the  harmonic  mean  2ab/(a  +  b)  of 
the  endpoints  of  Z^  ”  [a,b],  and  a  bound  for  its  percentage  error  is  then  100 
times  the  relative  width  (b  -  a)/(a  +  b)  of  ZL  (30).  The  maximum  percentage 
error  of  these  approximations  in  the  cases  in  which  existence  was  proved  by  the 
interval  iteration  is  shown  in  Table  8.2.  In  all  of  these  cases,  accuracy  to 
better  than  10%  was  achieved,  and  thus  the  results  are  guaranteed  to  be  of 
accuracy  which  is  sufficient  for  many  practical  purposes,  or  what  is  sometimes 
called  "engineering  accuracy".  This  accuracy  was  not  increased  by  requiring  the 
interval  iteration  to  converge  to  more  decimal  places. 

The  calculations  were  carried  out  on  a  Z80~based  microcomputer,  using  the 
Paacal-SC  compiler  of  Kulisch  and  Wippermann  [16],  which  provides  accurately 
rounded  floating-point  and  interval  arithmetic.  The  results  clearly  indicate 
that  by  going  to  "higher  resolution"  (more  subintervals  of  [0,1]),  the  range  of 
X  for  which  existence  of  the  solution  of  (8.1)  can  be  proved  can  be  extended, 
and  the  accuracy  of  the  numerical  approximations  obtained  can  be  increased. 

Computed  results  are  given  in  Appendix  I  for  X  “  0.92,  where  it  is  shown 
that  the  initial  interval  XQ  «  [1,1.75]  does  not  contain  a  solution  of  the 
integral  equation,  the  interval  iteration  converges  for  the  initial  interval  XQ 
-  [1,2.00],  but  the  existence  test  is  inconclusive,  while  XQ  “  [1,2.25]  is 
proved  to  contain  a  solution  of  the  integral  equation  which  is  bounded  by  the 
converged  values  of  the  interval  iteration.  Appendix  II  gives  the  source  code 
of  the  Pascal-SC  program  CHANSEK5  with  which  the  calculations  were  performed. 

The  microcomputers  used  were  a  Zilog  MCZ-1  with  the  RIO*  2.06  operating  system 
and  an  IMS  8000  SX  with  the  CP/M"  2.24  operating  system.  Both  systems  gave 
identical  results. 
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Table  8.1  (Continued) 


Initial 

Interval 


Lambda 


Existence  Convergence  Nonexistence 


(1,3.25] 


(1,3.50] 


[1,3.75] 


Existence  established  in  4  iterations  using  20  subdivisions  of 


•  undefined  for  the  initial  interval. 


Table  8.2 


Initial  Interval  I  Lambda  Iterations  Maximum  Percentage  Error 


References 


1.  A.  Abian.  The  Theory  of  Sets  and  Transfinite  Arithmetic.  W.  B.  Saunders, 

Philadelphia,  1965. 

2.  G.  Alefeld.  Intervallanalytische  Methoden  bei  nichtlinearen  Gleichungen. 

Jahrb.  Oberblicke  Math.  1979  (1979),  63-78. 

3.  G.  Alefeld  and  J.  Herzberger.  Einftihrung  in  die  Intervallrechnung. 

Bibliographisches  Institut,  Mannheim/Wien/Zttrich,  1974. 

4.  G.  Alefeld  and  J.  Herzberger.  Introduction  to  Interval  Computations 

(Jon  Rokne,  translator).  Academic  Press,  New  York  1983. 

5.  G.  Birkhoff.  Lattice  Theory.  AMS  Colloquium  Publications,  vol.  25,  rev. 

ed. ,  American  Mathematical  Society,  New  York,  1948. 

6.  0.  Caprani  and  K.  Madsen.  Contraction  mappings  in  interval  analysis. 

BIT  15  (1975),  362-366. 

7.  O.  Caprani  and  K.  Madsen.  Mean  value  forms  in  interval  analysis. 

Computing  25  (1980),  147-154. 

8.  O.  Caprani  and  K.  Madsen.  Interval  methods  for  the  solution  of  integral 

equations.  Report  No.  NI-80-05,  Numerisk  Institut,  Danmarks  Tekniske 

H*(jskole,  Lyngby,  1980. 

9.  O.  Caprani,  K.  Madsen,  and  L.  B.  Rail.  Integration  of  interval  functions. 

SIAM  J.  Math.  Anal.  12  (1981),  no.  3,  321-341. 

10.  S,  Chandrasekhar.  Radiative  Transfer.  Dover,  New  York,  1960. 

12.  J.  Garloff.  Interval  Mathematics:  A  Bibliography.  Freiburger  Intervall- 
Berichte  82/3  (1982),  169pp. 

12.  J.  L.  Kelley.  General  Topology.  Van  Nostrand,  New  York,  1955. 

13.  R.  Krawczyk.  Newton-Algorithmen  zur  Bestimmung  von  Nullstellen  mit 

Fehlerschranken.  Computing,  4  (1969),  187-201. 

14.  R.  Krawczyk.  Interval  extensions  and  interval  iterations.  Computing  24 

(1980),  119-129. 

15.  U.  Kulisch  and  W.  L.  Miranker.  Computer  Arithmetic  in  Theory  and 

Practice.  Academic  Press,  New  York,  1981. 

16.  U.  Kulisch  and  H.-H.  Wippermann.  Pascal-SC:  Pascal  for  Scientific 

Computation,  Institute  for  Applied  Mathematics,  University  of 
Karlsruhe,  and  Department  of  Computer  Science,  University  of 
Kaiserslautern,  1980. 

17.  R.  E.  Moore.  Interval  Analysis.  Prentice-Hall,  Englewood  Cliffs,  N.  J., 

1966. 


-23- 


.1 


18.  R.  E.  Moore.  Methods  and  Applications  of  Interval  Analysis.  SIAM  Studies 

in  Applied  Mathematics,  2,  Society  for  Industrial  and  Applied 
Mathematics,  Philadelphia,  1979. 

19.  R.  E.  Moore  and  S.  T.  Jones.  Safe  starting  regions  for  iterative  methods. 

SIAM  J.  Numer.  Anal.  14  (1975),  1051-1065. 

20.  K.  Nickel.  On  the  Newton  method  in  interval  analysis.  MRC  Technical 

Summary  Report  #1136,  University  of  Wisconsin-Madison,  1971. 

21.  K.  Nickel.  Stability  and  convergence  of  monotonic  algorithms.  J.  Math. 

Anal.  Appl.  54  (1976),  157-172. 

22.  L.  B.  Rail.  Numerical  integration  and  the  solution  of  integral  equations 

by  the  use  of  Riemann  sums,  SIAM  Review  7  (1965),  no.  1,  55-64. 

23.  L.  B.  Rail.  A  note  on  the  convergence  of  Newton's  method.  SIAM  J. 

Numer.  Anal.  11  (1974),  no.  1,  34-36. 

24.  L.  B.  Rail.  Computational  Solution  of  Nonlinear  Operator  Equations. 

Wiley,  New  York,  1969.  Reprinted  by  Krieger,  Huntington,  N.  Y.,  1979. 

25.  L.  B.  Rail.  A  comparison  of  existence  theorems  of  Kantorovich  and 

Moore,  SIAM  J.  Numer.  Anal.  17  (1980),  no.  1,  148-161. 

26.  L.  B.  Rail.  Integration  of  interval  functions  II.  The  finite  case. 

SIAM  J.  Math.  Anal.  13  (1982),  no.  4,  690-697. 

27.  L.  B.  Rail.  Solution  of  finite  systems  of  equations  by  interval 

iteration.  BIT  22  (1982),  233-251. 

28.  L.  B.  Rail.  A  theory  of  interval  iteration.  Proc.  Amer.  Math.  Soc.  86 

(1982),  no.  4,  625-631. 

29.  L.  B.  Rail.  Mean  value  and  Taylor  forms  in  interval  analysis.  SIAM  J. 

Math.  Anal.  14  (1983),  no.  2,  223-238. 

30.  L.  B.  Rail.  Representations  of  intervals  and  optimal  error  bounds.  Math. 

Comp.  41  (1983),  no.  163,  219-227. 

31.  L.  B.  Rail.  Application  of  interval  integration  to  the  solution  of 

integral  equations.  J.  Integral  Equations  5  (1983)  (to  appear). 
Preprint:  MRC  Technical  Summary  Report  #2128,  University  of  Wisconsin- 
Madison,  1980. 

32.  J.  T.  Schwartz.  Nonlinear  Functional  Analysis.  Gordon  and  Breach,  New 

York,  1969. 

33.  w.  Sierpifiski.  General  Topology  (C.  Cecilia  Krieger,  translator).  Univ. 

of  Toronto  Press,  Toronto,  1952. 

34.  P.  Wisskirchen.  Vergleich  iatervallarithmetischer  Iterationsverfahren. 

Computing  14  (1975),  45-49. 


-24- 


APPENDIX  I 


Computed  Results  for  X  =  0.92 


OXQPC  CHANSEK5 

ENTER  INITIAL  INTERVAL 

*[1,1.75] 

ENTER  LAMBDA 
*0.92 

EMPTY  INTERSECTION  AT  ITERATION  NUMBER  3 
FOR  LAMBDA  -  0.92, 

INTEGRAL  EQUATION  HAS  NO  SOLUTION  IN  THE  INTERVAL  HO  -  [  1.000,  1.750] 


H  [  1] 

-  [ 

1.00000000000, 

1.19034842076] 

2] 

-  [ 

1.16381274032, 

1.32334334049] 

H  [  3] 

-  [ 

1.27655952657, 

1.43596418062] 

H  [  4] 

-  t 

1.37089212268, 

1 .53512110298  ] 

H[  5] 

-  [ 

1.45333391818, 

1.62409194952 ] 

H  [  6] 

-  [ 

1.52700493008, 

1.70485377468] 

!  H[  7] 

-  [ 

1.59374951728, 

1.75000000000] 

1  H  [  8] 

-  [ 

1.65476948954, 

1.75000000000] 

H[  9] 

-  [ 

1.71089034579, 

1.75000000000] 

H  [  10  ] 

-  [ 

1.76270273459, 

1.75000000000] 

♦♦IMPROPER  INTERVAL** 

ENTER 

*1 

ENTER 

*[1,2 

ENTER 

"I"  FOR  NEW  INTERVAL,  "L"  FOR  NEW 

INITIAL  INTERVAL 

001 

LAMBDA 

LAMBDA,  "Q"  TO  QUIT 

*0.92 

FOR  LAMBDA  -  0.92, 

INTERVAL  ITERATION  CONVERGED  AT  ITERATION  NUMBER  10  TO: 


H  [ 

1] 

m 

l 

1.00000, 

1.19869] 

H  [ 

2] 

m 

1 

1.16670, 

1.34164] 

H[ 

3] 

m 

1 

1.28205, 

1.46512] 

H  [ 

4] 

m 

[ 

1.37845, 

1.57570] 

H  [ 

5] 

m 

t 

1.46232, 

1.67643] 

H[ 

6] 

m 

[ 

1.53677, 

1.  76911  ] 

H  [ 

7] 

m 

[ 

1.60369, 

1.85498] 

H  [ 

8] 

m 

[ 

1.66437, 

1.93493] 

H  [ 

9] 

m 

[ 

1.71977, 

2.00000] 

H  [ 

10] 

m 

[ 

1.77063, 

2.00000] 

♦♦EXISTENCE  OF  SOLUTION  NOT  GUARANTEED** 

ENTER  "I"  FOR  NEW  INTERVAL,  "L"  FOR  NEW  LAMBDA,  "Q"  TO  QUIT 
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*1 

ENTER  INITIAL  INTERVAL 

*[1,2.251 

ENTER  LAMBDA 

*0.92 

FOR  LAMBDA  -  0.92, 

HYPOTHESES  OF  SCHAUDER  THEOREM  VERIFIED  AT  ITERATION  NUMBER  2 
SOLUTION  OF  INTEGRAL  EQUATION  IS  BOUNDED  BY: 

H[  1]  -  [  1.00000,  1.22918] 

H[  2]  -  [  1.15545,  1.39469} 

H(  31  -  [  1.26068,  1.53896] 

H{  4]  -  [  1.34893,  1.66824] 

H[  51  -  [  1.42696,  1. 78559] 

Hi  6 ]  -  I  1.49798,  1.89285] 

H[  7]  -  I  1.56386,  1.99104] 

H [  8]  -  [  1.62568,  2.08055] 

H[  9]  -  [  1.68405,  2.16095] 

H[ 10 ]  -  [  1.73910,  2.23799] 

DO  YOU  WANT  TO  ITERATE  (Y/N)7 
*Y 

FOR  LAMBDA  -  0.92, 

INTERVAL  ITERATION  CONVERGED  AT  ITERATION  NUMBER  18  TO: 


H  [ 

11 

m 

[ 

1.00000, 

1.20069] 

H  [ 

2] 

m 

I 

1 .16670, 

1.34595] 

H  [ 

3] 

m 

[ 

1.28206, 

1.47194] 

Hi 

4] 

m 

1 

1.37845, 

1.58519] 

H  [ 

5] 

m 

( 

1.46232 , 

1.68868] 

H[ 

6] 

m 

I 

1.53677, 

1.78419] 

H  [ 

7] 

m 

l 

1.60369, 

1.87291] 

H  [ 

8] 

m 

[ 

1.66437, 

1.95573] 

H[ 

9] 

m 

[ 

1.71977, 

2.03334] 

H  [ 

10] 

m 

[ 

1.77063, 

2.10628] 

★♦EXISTENCE  OF  SOLUTION  GUARANTEED** 

ENTER  "I"  FOR  NEW  INTERVAL,  "L"  FOR  NEW  LAMBDA,  "Q"  TO  QUIT 
*Q 

KL/P-STOP 
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APPENDIX  II 


Source  Code  for  the  Pascal-SC  Program  CHANSEK5 


C> 

PROGRAM  CHANSEK5(INPUT, OUTPUT) ; 

CONST  DIM  -  10; 

TYPE  DIMTYPE  -  1 . .DIM; 

INTERVAL  -  RECORD  INF, SUP:  REAL  END; 

IVECTOR  -  ARRAY [ DIMTYPE ] OF  INTERVAL; 

RVECTOR  -  ARRAY[DIMTYPE30F  REAL; 

RMATRIX  -  ARRAY (DIMTYPE ] OF  RVECTOR; 

VAR  C:  CHAR; ALTER, ITER:  INTEGER ; HTEST , SUM , LAMBDA :  REAL;  I, 
HINF.HSUP:  RVECTOR ; LO , H I :  RMATRIX ; HO , COEF :  INTERVAL; 
GO, EXIST, EMPTY:  BOOLEAN; 

FUNCTION  INTPT  (  RA:  REAL  )  :  INTERVAL; 

EXTERNAL  41; 

OPERATOR  /  (  A , B :  INTERVAL  )  RES:  INTERVAL; 

EXTERNAL  85; 

FUNCTION  ILN  (  Y:  INTERVAL  )  :  INTERVAL; 

EXTERNAL  107; 

PROCEDURE  IREAD  (  VAR  F:  TEXT;  VAR  A:  INTERVAL  ); 

EXTERNAL  92; 


BEGIN  (*MAIN  PROGRAM*) 

FOR  J:-l  TO  DIM  DO  (*Generate  coefficient  matrices*) 
BEGIN 

COEF:-ILN( INTPT ( J+l ) /INTPT (J ) ) ;  (*Flrst  rows*) 

HI[1 , J] : -COEF. SUP ;LO[ 1 ,J] : -COEF.  INF; 

END;  (*First  rows*) 

FOR  I:-2  TO  DIM  DO 

BEGIN  (*Rows  2,. ...DIM*) 

FOR  J:-l  TO  (DIM-1)  DO 
BEGIN 

HI[I,J]:-HI[I-l,J+l];LO[I,J] : -LO[ I- 1 , J+l ] ; 

END; 

COEF: -ILN( INTPT ( DIM+I)/ INTPT< DIM+I-1) ) ; 
HI[I,DIM]:-COEF.SUP;LO[I,DIM]:«COEF.INF; 

END;  (*Matrix  generation*) 

C: • ' I '  ; WHILE  C  -  'I*  DO  (*Restart  with  new  initial  interval*) 
BEGIN 

WRITELN( ' ENTER  INITIAL  INTERVAL’); 

IREAD(INPUT,HO); 

C :  - '  L  ’ ; WHILE  C  -  'L'  DO  (*Restart  with  new  LAMBDA*) 

BEGIN 

WRITELN( ’ENTER  LAMBDA' );READ(LAMBDA); 

FOR  I : -1  TO  DIM  DO 
BEGIN 

HINF [ I J :-HO.INF;HSUP(I] : -HO. SUP 
END; 
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:  DIMTYPE; 


I 


HTFST :  -HO  .  SUP } 

GO: -TRUE } EMPTY  :-FALSE;EXIST:-FALSE;ITER:-l;VHILE  GO  DO 
BEGIN  ( *GO  Loop* ) 

(♦Update  HSUP [ 1  ]*  ) 

ALTER : “1 ; 

SUM :-l-> (LAMBDA/ >(2*DIM)*>SCALP( HI [ 1] ,HSUP,1) )} 

SUM  s “1 / >SUM ; 

IF  SUM  >  (HSUP[ 1 1-5.0E-06)  THEN  ALTER : -ALTER+1 } 

IF  SUM  <  HSUP [ 1 ]  THEN  HSUP ( 1 ] s -SUM} 

(♦Update  HINF [ I j  ,  HSUP [ I  ]  ,  I  -  2  .  .DIM*) 

FOR  Is-2  TO  DIM  DO 
BEGIN 

SUM : -1 - < (LAMBDA*  <(I-1)/<(2*DIM))*<SCALP( LO [I-l]fHINF,-I) ; 

SUM  s -1/ <SUM ; 

IF  SUM  <  (HINF[ I ]  +  5 .0E-06)  THEN  ALTER : -ALTER+1 ; 

IF  SUM  >  HINF[I]  THEN  HINF [ I ] : -SUM } 

SUM :  -1  -  > ( LAMBDA*  > I / > ( 2*DIM ) )*>SCALP( HI [ I ] ,  HSUP , 1) } 

SUMs-l/>SUM; 

IF  SUM  >  (HSUP[  H-5.0E-06)  THEN  ALTER : -ALTER+1 ; 

IF  SUM  <  HSUP [ I ]  THEN  HSUP l I ]: -SUM ; 

IF  HSUP [ I ]  <  H INF [ I ]  THEN  EMPTY : -TRUE } 

END;  (*Upda t e  of  HINF ( I ] , HSUP [ I  ]  for  I  -  2.. DIM*) 

IF  EMPTY  THEN 

BEGIN  (*Empty  Intersection*) 

WRITELN( 'EMPTY  INTERSECTION  AT  ITERATION  NUMBER  ' , ITER: 3) ; 
HRITELN ; WRITELN ( 'FOR  LAMBDA  -  ' , LAMBDA : 5 i 2 , ' , ' ) ; 

WRITELN( ' INTEGRAL  EQUATION  HAS  NO  SOLUTION 
IN  THE  INTERVAL  HO  -  [ ' , HO . INF i 6 s 3 , ' , ' , HO . SUP : 6 : 3 , ' ] ’ ) } 

HRITELN ; FOR  Is-1  TO  DIM  DO 
BEGIN  (*Improper  Interval  output*) 

WRITE( ' H I ' ,1:2, ' ]  -  [ ' , HINF ( I ] : 1 4 : 1 1 : - 1 , ' , ' ,HSUP[ I] :14: 11: 1 , '  ]  ’ )  ; 
IF  HSUP [ I ]  <  HINF t I]  THEN  WRITE('  **IMPROPER  INTERVAL**')} 
HRITELN} 

END}  (*Inproper  interval  output*) 

GO: -FALSE; 

END;  (*Empty  intersection*) 

IF  HSUP [ DIM j  <  HTEST  THEN 
BEGIN  (*Existence*) 

EXIST : -TRUE ; 

HRITELN ; HRITELN( 'FOR  LAMBDA  -  ' , LAMBDA : 5 : 2 , ' , ' ) ; 

HRITELN( 'HYPOTHESES  OF  SCHAUDER  THEOREM  VERIFIED  AT 
ITERATION  NUMBER  ',ITER:3); 

HRITELN( 'SOLUTION  OF  INTEGRAL  EQUATION  IS  BOUNDED  BY :'); HRITELN ; 
FOR  I : -1  TO  DIM  DO  (*Output  of  results*) 

BEGIN 

HRITELN( ' H [ ' ,1:2 , ' ]  -  [ '  , HINF [ I ] : 8 : 5  :  - 1 , '  , '  ,  HSUP [ I]  :  8  :  5  : 1 ,  '  ]  '  ) 
END}  (*Output  of  results*) 


WRITELN ; WRITELN ( 'DO  YOU  WANT  TO  ITERATE  ( Y/N ) ? ’ ) { READ ( C , C )  ; 

IF  C  <>  'Y*  THEN  GO: “FALSE  ELSE  HTESTs-O; 

END;  (*Existence*) 

IF  ALTER  -  20  THEN 
BEGIN  ( ^Convergence* ) 

GO: “FALSE ; 

WRITELN; WRITELN( 'FOR  LAMBDA  •  '  , LAMBDA : 5 : 2 , * , ' )  ; 

WRITELN( ' INTERVAL  ITERATION  CONVERGED  AT  ITERATION 
NUMBER  ' , ITER :  3 , '  TO: ’ )  ; 

WRITELN; 

FOR  I:“l  TO  DIM  DO  (*Output  of  results*) 

BEGIN 

WRITELN< • H[ * , 1 : 2 , * ]  -  f * ,HINF[I]:8:5:-1, * , ' , HSUP[1 ] : 8 : 5 : 1 , ’ ) ' )  ; 

END;  (*0utput  of  results*) 

WRITELN; 

IF  EXIST  THEN  WRITELN( ' **EXI$TENCE  OF  SOLUTION  GUARANTEED**') 

ELSE  WRITELN( '**EXISTENCE  OF  SOLUTION  NOT  GUARANTEED**'); 

END;  ( *CONVERGENCE* ) 

ITER : -ITER+I ; 

END;  (*GO  Loop*) 

WRITELN; 

WRITELN( 'ENTER  "I"  FOR  NEW  INTERVAL,  "L"  FOR  NEW  LAMBDA,  "Q"  TO  QUIT') 
READ(C,C) ; 

END;  (*Rester t  with  new  LAMBDA*) 

END;  (*Restart  with  new  Initial  Interval*) 

END.  ( *MAIN  PROGRAM*) 
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